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Summary. — Gigantic neutrino telescopes are primarily designed to search for very 
high energy neutrino radiation from the cosmos. Neutrinos travel unhindered over 
cosmological distances and therefore carry unique undistorted information about its 
production sites: the most powerful accelerators of hadrons in nature. 
In these lectures, we present the relevant physics motivations and their specifics. We 
review methodological aspects of neutrino telescopes: the experimental technique, 
some of the faced problems and the capabilities in terms of discovery potential, ef- 
fective area, isolation of a signal from atmospheric backgrounds, etc. Instruments 
and their operation in various media are described. We also mention the instrumen- 
tal birth and provide an outlook of the detection technique toward very low and 
ultra-high energies. 
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1. Introduction 

Following the discovery of cosmic rays (CR) in 1912, the field of elementary parti- 
cle physics was founded upon the identification of a number of related particles such as 
positrons, muons, pions and strange particles. One hundred years later, however, the 
origins of the CR themselves remain obscure. The prevalent bottom-up paradigm posits 
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their origin within cosmic accelerators, e.g. pulsar wind nebulae, X-ray binaries, super- 
nova remnants in our galaxy or active galactic nuclei and gamma ray bursts, via the so 
called " Fermi acceleration" mechanism or alternatively in large scale electric fields present 
in pulsar environments. The interaction of the accelerated particles with the surrounding 
radiation fields or dense gas clouds in the neighborhood produces a secondary neutrino 
flux. This leads to the experimental possibility of conducting astronomy through the 
means of high energy (HE) (> O(10 GeV)) neutrino messengers, perhaps providing the 
key to unlocking the mysteries of the CR origin. 

The nascent field of neutrino astronomy, at somewhat lower energies, has grown from 
the first successful observation of extra-terrestrial (ET) neutrinos from a supernova in 
1987 and the sun, bearing witness to the death of a star and precipitating the monumental 
achievement of establishing a neutrino mass. 

The first neutrino telescopes in operation, NT200+ at Lake Baikal and AMANDA at 
the South Pole, attempted to open a new window on the Universe from the perspective 
of the HE neutrino, successfully detected HE neutrinos of atmospheric origin and have 
motivated the construction of larger detectors, the IccCube detector at the South Pole 
instrumenting a cubic kilometer of ice and Antares in the Mediterranean Sea, both in 
operation. As yet, however, no high energy extra-terrestrial signal has been detected, 
but phcnomcnological arguments suggest that signal detection should be possible with 
the completed IceCube array after a few years of data taking. 

These detectors are not limited to the question of the CR origin and can be utilized 
to probe questions of fundamental particle physics, on the nature of the neutrino or on 
the validity of Nature's symmetries by means of the atmospheric HE neutrino beam, 
and cosmological scenarios. At the interface of these three pillars, fundamental particle 
physics, astrophysics and cosmology, HE neutrino telescopes may also be instrument in 
unveiling the nature of dark matter. This is discussed in Sections 2-6. 

Section 7 introduces the detection technique and several methodological aspects: sig- 
natures, notions of effective areas and discovery potential, analytical approaches, event 
reconstruction and analysis systematics. Section 8 presents historical elements in the 
experimental development of the field and its current status. 

The choice of topics described in Sections 9-11 is intended to cover a broad spectrum 
of the methodology. Section 9 shows that high energy neutrino telescopes can be utilized 
as supernova detection instruments, i.e. they can probe a neutrino signal at energies many 
orders of magnitude below their original concept. A back of the envelope calculation on 
their potential is presented. Section 10 describes the generic procedure for point source 
searches and introduces the important notions of upper limit, sensitivity and discovery 
flux. The results reported in recent years denote large methodological "breakthroughs" 
in the exploitation of these instruments with a sophistication of the analysis methods. 
Section 11 concerns the ultra high energy detection methodology with the presentation 
of the principle underlying alternative signatures. 

In these lectures, we do not emphasize actual experimental results, which can be 
found in the mentioned literature. 

2. Cosmic Accelerators and neutrino astronomy 

The generic picture of the cosmic ray origin appears quite satisfactory up to w 10 20 eV 
to explain the CR energy spectrum and composition changes: CRs are accelerated in 
the vicinity of galactic and extra-galactic accelerators through efficient first order Fermi 
mechanism [1]: particles, bouncing back and forth a strong shock propagating through 
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interstellar plasma, gain energy AE oc E at each crossing of the wave-front with a small 
probability of escape of the region. This mechanism trivially predicts a power law spectral 
shape 

(1) d$/d£ ~ E-i 

for the accelerated particles. Moreover, the spectral index takes a universal 7 = 2 for 
strong (supersonic) shocks with speed /3c, due to the specific momentum imprint to the 
ISM from the shock crossing [2, 3]. Roughly, nuclei with atomic charge Z will gain energy 
as long as they remain confined in the acceleration region of size L and magnetic field B 
{i.e. gyroradius smaller than L) and can therefore attain a maximum energy 

(2) £ max ~ ZpBL. 

Several outstanding features characterize the CR spectrum illustrated Fig. 1: 

- A first break in the power law with a steepening of the CR spectrum occurs at a 
few PeV (knee). This break, which may sign the limitations of the confinement of 
the galaxy or galactic accelerators [4, 5], is accompanied by a gradual change of 
composition from light to heavy nuclei. 

- A second break occurs at a few EeV (ankle) and is understood as a cross-over 
of the fading galactic component and a harder extra-galactic component. This 
interpretation is further supported by a change of composition from heavy to light 
nuclei and by the fact that the galaxy cannot confine nuclei at these energies. CR's 
with energy exceeding that of the ankle are called ultra high energy (UHE) CR's. 

- At ultimate energies, around w 1 o 20 ~ 21 eV per nucleon, protons interacting with the 
cosmic microwave background (CMB) undergo an energy damping due to photopion 
production over cosmologically short distances known as the GZK cutoff [6]. See 
also Section refsct:gzk 

Together with CR electrons, CR nuclei have provided valuable information in the 
study of prevailing acceleration processes and the morphology of the galaxy. However, 
the nature of cosmic accelerators is not revealed easily because CRs are poor astronomical 
messengers: deflected in magnetic fields, they do not point back to their sources. 

The occurrence of electronic acceleration (by the first order Fermi acceleration mech- 
anism) up to about 100 TeV was demonstrated indirectly in supernova remnants (SNR) 
by means of conventional astronomy observing gamma radiation extending up to tens 
of TeV: the accelerated electrons emit non thermal Bremsstrahlung radiation and syn- 
chrotron radiation, which in turn, together with other surrounding radiations become the 
targets for the inverse Compton process. The HE photons can then further produce elec- 
tron pairs by interacting with the surrounding radiation fields. This dynamics produces 
complex spectral energy distributions (SED), typically with a double humped structure. 

Following an old argument by Baade and Zwicky [7] the energy balance suggests that 
SNR's could account for the bulk of the galactic CR's (GCR) as well: given a CR density 
Pcr ~ 1 eV/cm 3 and a CR galactic confinement time of about t osc w 10 7 yr in a galactic 
volume Vj»ai w 10 67 cm 3 , the CR luminosity is 

(3) £cr = PcR^gai/T-csc ~ 10 41 erg/s. 
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With one supernova releasing 10 51 erg every 30 yr (10 9 s), £cr can be sustained provided 
a plausible conversion efficiency of about 10% of the SNR energy into CR kinetic energy. 

However, this postulate has not been unambiguously established through multi-wave- 
length (MWL) studies, because of the difficulty in disentangling the electronic component 
from the 7r°-decay component signing hadronic acceleration (from processes pp, P7 — > 
7r — > 7, u). The observation of HE neutrinos would provide an incontrovertible proof of 
hadronic acceleration. 



Neutrinos as astronomical messengers 
present further advantages, as they may es- 
cape optically thick sources and provide an 
unrestricted access to the distant HE sky. Dis- 
tant photons on the contrary interact with 
the extragalactic background light, produc- 
ing electronic pairs and nearby blazar spectra 
are distorted if not absorbed above TeV ener- 
gies^). However, neutrinos, with their small 
interaction cross sections are difficult to detect 
but searches for ET neutrinos are greatly alle- 
viated by the combination of an increase of the 
cross section with growing energy and rather 
hard expected flux spectra d$ v /dE v cx E~ 2 
as suggested by the first order Fermi accel- 
eration mechanism. This is the reason for 
building large HE neutrino telescopes, in par- 
ticular, IceCube [8] and Antares [9], already 
mentioned, currently the most sensitive in- 
struments in operation with complementary 
fields of views. The discovery of a positive 

signal through the exploration of the neutrino Fi S- L " The all-particle cosmic ray spectrum 
, u i i i • -c 4- u i and its main features, 

sky would clearly represent a significant break- 
through and begin to resolve a hundred year 
old mystery. This is the key ambition of neutrino astronomy. 




We note that conventional astronomy is continually revealing an ever richer sky at 
GeV and TeV energies, powered by astronomical objects. Among the favorite GCR ac- 
celerator candidates, besides Pulsar Wind Nebulae and Shell-like SNRs, which exhibit 
rather steady behavior, we find other objects, notably associated to binary systems, 
which exhibit variable gamma-ray emission [10, 11, 12, 13, 14, 15, 16]. Prime candidates 
to explain cosmic rays at the highest energies, beyond the ankle, are Active galactic 
nuclei (AGN), often displaying variability [17] and cataclysmic gamma ray burst (GRB) 
phenomena, which are quite common in the universe. The transient behaviors of some 
of these neutrino source candidates may facilitate the rejection of the atmospheric back- 
ground and is therefore interesting from an experimental standpoint; this is particularly 
true for GRB's with analyses essentially background free because of the localized emis- 



( ) The situation is even more dramatic for PeV gammas interacting with the Cosmic Microwave 
Background (CMB): if emitted in the galactic center, they would not reach us. 
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sion of a short burst of neutrinos, assuming the connection between gamma-rays and 
neutrinos discussed below. 

2'1. UHE CR cosmic ray connection. - The energy density of the extra-galactic com- 
ponent extracted from the green line on Fig. 1 is Puhecr ~ 3 • 10~ 19 erg/cm 3 . Therefore, 
the required power for a source population to generate puhecr over a Hubble time of 
*h = 10 10 yr (per unit Mpc 3 ) is 

(4) £uhecr « 10 45 erg Mpc~ 3 yr" 1 . 

From the perspective of astrophysics, the coincidence between this number and the elec- 
tromagnetic output from AGN (2 • 10 44 erg/s per AGN) and/or GRB (2 • 10 52 erg per 
GRB) source populations, which fits the generic hypothesis of a transparent source, is the 
reason why these objects have come to be among the favored extra-galactic candidates 
for UHE hadronic acceleration. 

In the transparent source, while accelerated protons remain trapped in the acceler- 
ation region and interact with the surrounding radiation field, the secondary neutrons 
and neutral decay products of secondary pions are no longer confined. This results in a 
similar energy injected into CR, 7 and v escaping the source and related spectra: the P7 
reaction proceeds through various modes, direct production, resonant (and multi-pions) , 



(5) P7 — > 7T + n, 7T°p 

and the neutrons, pions and muons further decay, 

(6) n — > 7r~p 

(7) 7r ± -> p^Vy, -> e ± u ll v ll v c 

(8) 7T° -> 77 



Near the photo-pion production threshold, assumed to dominate, pions take away about 
TO^/TOp energy fraction, while asymptotically pions take 50%. Charged pions decay 
kinematics is distributing energy roughly a quarter to each final state particles. Even- 
tually, the same amount of energy goes into gamma and neutrinos. For pp interaction, 
the picture brings to similar conclusion. Therefore, assuming a E~ 2 neutrino spectrum, 
(/ E V <3& V /&E V = cp UH ecr/47t), we have, 

~ 5 ' 10- 8 GeVcm 2 s" 1 sr" 1 . 

aE v 

The diffuse neutrino flux from this generic transparent CR source is referred to as the 
Waxman-Bahcall (WB) upper bound [18], once normalized on the UHECR component, 
and indicates the necessity to consider neutrino telescopes at the cubic kilometer scale 
(see Section 7). This upper bound can be relaxed by assuming some optical opacity 
of the source, while maintaining consistency within experimental constraints (CR data 
and diffuse gamma-ray background) . This is known as the Mannheim-Protheroe-Rachen 
(MPR) upper bound [19]. 

We discuss now the blazar as an example of hadronic acceleration to illustrate further 
the relation between gamma ray and neutrino emission. 



(9) 
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2'2. AGN blazar accelerator. - Blazars, an extreme manifestation of AGNs when 
the jets are pointing closely toward the observer are strong candidate sites for EG CR 
acceleration and for subsequent neutrino emission within hadronic models of activity 
given the compelling arguments for non thermal emission. Blazars come in two flavors, 
Flat Spectrum Radio Quasars (FSRQ) which shows broad emission line regions and low 
to high frequency BL Lac (LBL — > HBL) characterized by the absence of broad emission 
lines and a smaller luminosities. A characteristic double humped spectrum is observed, 
which peaks in the infrared & MeV bands for FSRQ and in the optical/UV - UV/X- 
ray & GeV-TeV bands for LBL and HBL. Many models of blazar activity predict the 
associated neutrino flux within reach of the current generation of instruments. They 
exhibit variability at all wavelengths and on various time scales (flares) which can be 
as short as one hour, constraining the size of the emission region to the submilli pc 
scale. A few blazars have shown remarkable and unusual behavior, emitting orphan 
flares, i.e. TeV outburst emissions with no counterpart at any other wavelengths [20]. 
This peculiar observation is hardly reconcilable with leptonic models of activity which 
predict a correlation between X-ray and TeV emission. Instead, it could sign the decay of 
mesons following p7 interaction of accelerated hadrons, which then result in the emission 
of neutrinos [21] (the constraints on the correlated emissions are looser in hadronic models 
of activity). 

There are two broad categories of models of blazar activity, leptonic and hadronic 
and both can explain the double humped spectrum. A realistic model of 7-ray emission 
is likely to be the synthesis of both, if these objects are to be at the origin of a significant 
fraction of the extra-galactic CR flux. 

Within the leptonic class of scenarios, protons are not accelerated at energies high 
enough for efficient P7 interaction and the radiation is dominated by relativistic elec- 
trons. Synchrotron Self Compton (SSC) and Synchrotron External Compton (SEC) 
model classes are often considered in the literature [22]: external radiation, e.g. from the 
accretion disk, or synchrotron radiation within the jet of accelerated electrons constitute 
the target for the inverse Compton process in SEC resp. SSC models. 

Within the hadronic class of scenarios, when a significant fraction of the jet kinetic 
power is converted into the acceleration of relativistic protons, two main model classes 
are distinguished based on P7 or pp interactions. 

In the Synchrotron Proton Blazar model (SPB) [23], protons are co-accelerated with 
electrons but the energy density of electrons is small compared to the energy density 
carried by accelerated protons. However, the synchrotron emission by electrons accounts 
for a large fraction of the X-ray radiation. The main TeV radiation process from pro- 
ton synchrotron radiation is an essential ingredient in the reproduction of the double 
humped structure of the SED. Multiple processes initiate a HE 7 radiation, interaction 
of protons with the X-ray target, P7 Bethe-Heitler pairs, p, tt and fi synchrotron losses 
and electrons from muon decay, inducing EM cascades, which redistribute the power to 
lower energies. Eventually, the emission region becomes optically thin and the radiation 
escapes. Contrary to the Proton Initiated Cascade model (PIC) [24], which predicts 
large neutrino fluxes but is unable to reproduce the complex SED structure, very large 
magnetic fields are necessary in the SPB model and the resulting neutrino flux can be 
much lower. Experimentally, this model is interesting if periods of enhanced activity are 
considered (e.g. Mrk 501 in 1997) as the hardening of the spectra during these periods 
does not only translate into a larger neutrino flux but also into evolving spectral features. 

Within the framework of pp scenarios [25], where accelerated protons interact with 
cold protons present in the jets, the observed double humped SED structure can be 
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reproduced as well. In addition, given the low energy threshold of the reaction, this type 
of models requires only relatively low magnetic fields. The interaction rate can be high, 
resulting mainly in pions, which on decay produce neutrinos, the MeV-TeV radiation 
and further synchrotron emission from EM cascades may account for the observed X-ray 
radiation. The bolometric luminosity is related to the neutral pion component. 

In the hadronic models of activity [26], the maximal flux of neutrinos allowed for pp 
and p7 interactions was presented. In the case of an isotropic injection of soft protons 
(Fermi accelerated) and subsequent pp so ft interactions, the neutrino spectral shape fol- 
lows closely the HE 7 spectral shape, as a result of source transparency to TeV radiation. 
Applied to the list of blazars characterized with Fermi [27], the Fig. 2 shows the discovery 
potential of IceCube, which was analytically calculated in [26], in a <I> 7 vs T 7 diagram. 
Notice that the flux (integrated above 100 MeV) necessary for a detection drastically 
changes with the spectral index. Consequently, a few BL Lac would be within reach of 
IceCube but none of the bright FSRQs (due to softer emission spectra) often considered 
as favored neutrino source candidates in the literature. 

Within this model, these results convey the 
importance of the necessity for a "hybrid" ap- 
proach from the methodological standpoint in 
search for a neutrino emission, i.e. to consider 
the 7-ray integrated flux in combination 
with the corresponding spectral index in order 
to optimally search for neutrinos. Moreover, 
assuming that an enhanced neutrino emission 
and a hardening of the neutrino spectrum may 
be expected from blazars in high states, this 
prompts an analysis methodology, which relies 
on the time-dependence of the spectral shape 
and normalization of the neutrino emission,Fig. 2. - From [26]. Integrated flux $ 7 w.r.t. 
for optimizing the discovery potential. the spectral index T 7 distribution of Fermi 

AGN [27] . Superimposed, the gray curves rep- 

r . 11-111 resent the discovery threshold of IceCube in 

In [26], p 7 and pp models with hard pro- the ppsoft model (gee text) 

ton injection spectra were considered. In this 
case, the bolometric luminosity provides con- 
straints on the upper neutrino flux level but 

the neutrino flux spectral shape, due to cascading processes of TeV photons, is not di- 
rectly related to the GeV band spectrum. Therefore, in the eventuality of HE neutrino 
detection, the location of the detected sources in the diagram would strongly constrain 
the source acceleration mechanisms, opening the captivating prospects of delineating be- 
tween hard (e.g. large scale electric fields at the base of the jet) and soft (first order 
Fermi acceleration) acceleration mechanisms. 

We have presented above the basic motivations of neutrino astronomy and provided 
several hints calling for a generic point source search methodology. The reader will 
refer to the lectures by G. Sigl, in this school, for a more detailed discussion on the 
phenomenology of high energy neutrino sources. We limit ourself to the methodological 
aspects in the following. 
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3. Atmospheric neutrinos 

Along with muons, neutrinos are produced by the interaction of CRs with the atmo- 
sphere. Atmospheric neutrinos constitute an irreducible background in searches of ET 
neutrinos, contrary to atmospheric muons. The atmospheric neutrino spectrum can be 
split into a steeply falling conventional component [28, 29, 30] and into a hard charmed 
or prompt component [31, 32, 33, 34]. 

The conventional atmospheric neutrino component is relatively well known (uncer- 
tainties remain on the primary CR flux and the subsequent air shower modeling the 
kaon production) and results from the disintegration of pions and kaons, which progres- 
sively reinteract at growing energies. This "beam dump" effect increments the spectral 
index by one unit (7 = 2.7 — » a = 3.7) ( 2 ) This softening is essential in making the 
case of HE neutrino telescopes: eventually, there will be an energy above which a flux of 
ET neutrinos, assumed harder, will start to dominate. Kaon-induced neutrinos become 
increasingly dominant with energy and contribute already about 75% of the neutrinos at 
TeV. 

In contrast, the feeble prompt component remains hard up to PeV energies: the fast 
decay of very short lived charmed hadrons and mesons (produced as pairs DD, A C D, etc.) 
prevents the "beam dump" to be effective up to energies well above PeV. Therefore, the 
prompt flux eventually surpasses the conventional flux above some critical energy, with a 
poorly estimated threshold around 0.1 — I PeV due to the large uncertainties affecting its 
normalization from largely unknown cross sections (O(0.0I — 0.1 mb) in the range 0.1 — 1 
PeV). 

The angular distribution of the conventional and the prompt components are quali- 
tatively different: below the critical energy, the absence of a beam dump for the short 
lived charmed mesons and hadrons results in an angular prompt neutrino distribution 
independent of cos 9, where 9 is the zenith angle. For conventional mesons, the reinter- 
action rate is sensitive to the interplay between the decay length and interaction length, 
a function of cos 9: inclined showers develop at first in comparatively low density atmo- 
sphere and therefore, the mesons preferentially decay compared to vertical showers up 
to higher energies, therefore resulting in an increase of the HE neutrino flux. The flux 
suppression will depend linearly on cos#, reflecting the slant depth behavior. 

Therefore, an additional discrimination difficulty in searches for an isotropic diffuse 
flux of very high energy (above PeV) neutrinos arises from the isotropic prompt neutrino 
flux component, even if an ET neutrino flux is expected to be slightly harder. 

3T. CR muon background. - The muon intensity decays exponentially with slant 
depth (the integrated column density along a path) penetration in the ground. At depths 
of about 14 kmwc, it turns into a constant: the muon flux is induced by neutrinos 
interacting in the neighborhood of the detector. Neutrino telescopes, projected or in 
operation, are located at depths between 1 — 5 kmwe, i.e. at depths where the CR- 
induced muon flux is respectively 10 7 — 10 4 times larger than the neutrino-induced muon 
flux. CR-induced muons constitute a reducible background to neutrino searches, in 
particular to atmospheric neutrino study, which depends on the detector quality, i.e. its 



( 2 ) In analogy to the spectral softening of electrons propagating in the galaxy, the linear increase 
of the pion lifetime in the expanding shower together with the linear increase of the energy loss 
per track length with energy results in an energy loss before decay in E 2 , therefore the muon 
spectral index is a = 7 + 1. 
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ability to reconstruct the incoming lepton direction. Thus, it is clear that burying the 
detector as deeply as possible helps to greatly alleviate this task. For instance, in order to 
isolate a sample of upward-moving atmospheric neutrinos with 90% purity, the required 
discrimination power between downward muons and neutrinos is estimated at the 10 7 
level at a 2 kmwe depth. 

3'2. Strong ET neutrino discrimination. - HE neutrinos events are conjointly pro- 
duced with muons in CR showers, opening the possibility of discrimination between 
downward-moving extra-terrestrial and atmospheric neutrino events: the latter being 
accompanied by a muon (bundle), in particular a muon from the same parent meson. 
In [38] , it was actually shown that downward- moving atmospheric neutrino could be effi- 
ciently rejected in an IceCube-like detector with external VETO layers for zenith angles 
up to 60° and neutrinos with energy in excess of 10 TeV. 

The atmospheric neutrino beam is interesting on its own as a a signal for particle 
physics study (e.g. for charm production cross sections, neutrino oscillations, tests of 
the Lorentz violation, etc.). In principle it enables to picture the tomography of the 
earth [35]. 

4. Neutrino oscillations 

Similarly to the quark sector, mass and weak eigenstates do not coincide in the neu- 
trino sector. The unitary neutrino mass matrix U a i relates mass to weak eigenstates 
according to \u a ) — J2i U a iWi). The evolution of a neutrino produced in a certain flavor 
state a is derived from the evolution of its mass eigenstates, i.e. governed by 

(10) \v a (x,t)) =J2e- i ^ t ~^U ai \u i ) =Y. U <» e ~ l{mV2E)tU >fi) 

i i,0 

in the limit 7> 1. The propagated state differ in general from the initial state and the 
probability of appearance of a flavor /3 at distance L is 

(11) P^u,(L) = \{u p \v a {L))\ 2 = S aP - 4]Tsin 2 {5m%L/AE) U m U Pl U aj U Pj . 

j>i 

(considering real MNS matrix elements is correct for oscillations calculations in absence 
of CP violation). Experimentally and in "good" approximation, {J7 e ,}j w {ci 2 ,si 2 ,0}, 
{Upi}i « ^{-Si2,ci2,l} and {U Ti } l w ^{«12, -C12, 1}, with c i2 = cos#i 2 and s 12 = 
sin #i2, where 6 12 = 9 &oi w 33°, assuming # 23 = #atm ~ 45° and 613 w 0°; <5m 21 = Sm^ ol ss 
8 • 10- 5 eV 2 and \Sm 2 31 \ ss \6m 2 32 \ = 6ml tm ss 2.4 • lO" 3 eV 2 (<5m 2 3 = m 2 - m 2 ). 

We leave it to the reader using eq. 11 to verify that the first dip in the neutrino oscil- 
lation probability — > is at about 25 GeV for a neutrino crossing the earth through 
its center. Using the atmospheric neutrino beam, this region could be explored [36, 37] 
provided a neutrino telescope with a sufficiently low energy threshold. 

4T. Oscillations over astronomical baselines. - Distant astrophysical candidates of 
high energy neutrinos feature Sm 2 L/E 3> 1. This implies that the neutrino flux measured 
with a finite energy resolution AE/E is averaged over many periods of the oscillation 
probability argument sin 2 (5m 2 L/AE) -> / dE' sin 2 (5m 2 L / AE')g a (E' - E) = 1/2, where 
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g a (E' — E) is a Gaussian distribution with width a = AE/E for instance (moreover, 
the acceleration region is usually extended, of size AL). One can therefore drop the 
dependences in Sm 2 , E and L for the calculation of the neutrino flux at Earth, 



(12) 



S £>>1 = <W - 2 J2 u ai U pi U aj U t 



03- 



3>i 



Assuming the special 6*23 and $13 values and starting with a neutrino mixture at the 
production site 



(13) 



= {1,2,0}, 



which follows from both processes P7 and pp (summing neutrinos and anti-neutrinos) , it 
is straightforward to obtain the mixture propagated to earth 



(14) 



{4v.,^ M ,^ T } = {1,1,1}: 



all three flavors are equally populated at earth, opening the favorable experimental per- 
spective of potential source detection with distinct signatures. Moreover, the presence 
of tau neutrinos potentially enables to look at up-going events in the very high energy 
regime, otherwise hindered. 

If we would have kept the distinction between neutrinos and anti-neutrinos, we would 
have noticed that pp and P7 leads to distinct mixtures. It is left as an exercise to 
conclude that the D e flux is relatively lower in case of P7, opening the experimental 
possibility of discriminating the dominant process at work in a cosmic accelerator via 
the Glashow resonance, provided it is detected at PeV energies (relative rates of track-like 
and cascade-like events are affected, see Section 7). 

The pattern of neutrino oscillation is affected in the presence of matter or assuming 
violation of the Lorentz invariance and briefly outline these topics below. 

4'2. Matter oscillations. - In the pres- 
ence of regular matter, the interaction of the 
electronic flavor differs from the other fla- 
vors. While in regular matter v c CC scatter- 
ing on electrons is possible, this is forbidden 
for the other flavors and leads to the so-called 
Mikheyev - Smirnov - Wolfenstein (MSW) ef- 
fect. In more formal terms, an additional po- 
tential appears in the Hamiltonian for the CC 
interaction and uniquely for the electronic fla- 
vor, so that mass eigenstates in vacuum are 

no longer mass eigenstates in matter. Work-Fig. 3. - -> i/ M transition probability 
ing out the solution found in standard text- of upward-moving neutrinos for different 13 
books, we notice that the presence of matter values and hierarch y assumptions after earth 
eventually enables the disentanglement of the crossm ^' 
neutrino hierarchy for non vanishing sin 2 (9 13 

values (in vacuum we only have access to the absolute value of the difference of mass 
squared, not its sign). 
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A neutrino telescope with a O(10 GeV) energy threshold may resolve the hierarchy, 
based on the statistical discrimination of v from D (CC interaction cross sections differ 
by a factor two)( 3 ), provided rather large values of sin 2 6*13 > 0.1, as the oscillation 
pattern of neutrinos in the range between 10 and 15 GeV crossing the earth's core is 
most strongly affected by matter effects [39]: matter effects enhance (or suppress) the 
oscillation probability of neutrinos or anti-neutrinos if the hierarchy is normal or resp. 
inverted (inverted resp. normal). Fig. 3 illustrates the — > transition probability of 
vertically up-going neutrinos for different sin 2 #13 values and hierarchy assumptions after 
earth crossing (i.e. reaching the detector). 

From an experimental point of view, resolving the neutrino hierarchy may turn out to 
be a difficult analysis dominated by systematic uncertainties: good energy resolution and 
detection efficiency down to £?„ = 5 GeV are required (i.e. the ability of reconstructing 
tracks about 30 m long with about 10 m resolution), as well as a Monte Carlo provid- 
ing relative flux systematic uncertainty down to of a few percents w.r.t. the incoming 
neutrino direction on a moderately large angular scale (the kinematics enables only a 
neutrino incoming angle resolution of (6v M ) ~ 10° at these energies). 

It is worth mentioning that, in principle, the possibility to "distinguish" event by 
event between and by assigning a tag probability exists. This, without measuring 
the charge of the induced muon, but instead based on the precise topology of the con- 
tained event: the measurement of both, the energies of the induced muon (through its 
track length) and the energy of the induced shower at the interaction vertex, grants ac- 
cess not only to the incident neutrino energy but to the interaction inelasticity y as well. 
The average value taken by y differs for neutrinos and anti-neutrinos. Such an approach, 
in contrast to the statistical v — v discrimination based on cross section difference, will 
greatly help deciphering the neutrino hierarchy. 

4'3. Lorentz invariance violation. - Many efforts to construct a quantum gravity [40, 
41, 42] theory imply the deformation or the breaking of the Lorentz symmetry and pre- 
dict non trivial modifications of space-time symmetries at the Planck scale, such as the 
existence of a new fundamental length scale lp\. The Standard Model Extension frame- 
work [43] provides an effective field-theoretical approach for the study of the violation of 
Lorentz invariance. Among the most promising experimental signatures [44, 45, 46, 47] 
within this phcnomcnological framework are neutrino oscillations as flavor changing sig- 
natures are amplified. 

The modification of the dispersion relation, retaining the validity of the energy- 
momentum conservation law, is a simple kinematic framework to introduce a violation 
of the Lorentz invariance: 

(15) E 2 = m 2 + (l + f l )p 2 , 

where i denotes the energy eigenstate and where the {fi}i\fi <C 1 depend on the species. 
At energies above 100 GeV (in order to neglect mass- induced oscillation), the oscillation 
probability, left as an exercice, is 

(16) p v ^ Vfl = 1 - sin 2 26 v sin 2 (SvEL/2), 



( ) Neutrino telescopes cannot distinguish between the negatively and positively charged in- 
duced muons. 
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where 9 V is the maximum attainable velocities (MAV) mixing angle to the flavor eigen- 
states and 5v the fractional MAV difference. We note that (1) oscillations also occur if the 
neutrino mass vanishes and (2) the EL dependence of the argument instead of the usual 
L/E in standard oscillations. The exotic term competes with the classical term and may 
dominate at high energy, leading to distortions of the angular and energy distributions 
relative to pure mass-induced oscillations. Distribution discrepancies beyond statistical 
and systematic errors would suggest a violation of the Lorentz invariance. Alternatively, 
this would help constrain the allowed parameter space (5v, 9 V , etc.). IceCube in its 
baseline configuration will accumulate a unique sample of O(10 6 ) atmospheric neutrinos 
with energies above 50 GeV, representing an increase in sensitivity for testing a violation 
by more than an order of magnitude on existing constraints [48, 49, 50], in the particular 
case of coinciding mass and asymptotic velocity eigenstates [51]. A violation could be 
probed up to a fractional velocity difference approaching 5v = 10~ 28 , the limiting factor 
being the shrinking size of the neutrino sample above a growing energy threshold. 

5. GZK neutrinos 

UHE protons may interact with low energy photons (cosmic microwave or infrared 
background) P7 — > 7rn, provided the CoM energy of the reaction satisfies s > (m p + m v ) 2 . 
This corresponds to an energy threshold of w 10 20 eV for a frontal collision and a mean 
CMB photon energy of (e) = 6.4 • 10~ 4 eV and results in the degradation of the nucleon 
energy and the appearance of the so-called GZK cutoff [6]. 

Following the photo-pion production reaction, neutral and charged pions and neutrons 
are produced, which further decay into gamma rays and neutrinos. 

The interaction length is approximately given by 



where a pl — 0.25 mb and n 7 = 411 cm~ 3 . 

Accounting more precisely for the exact shape of the photo-pion production cross 
section a pl (E p , e 7 , cos 9) with direct, multi-pion, diffractive and resonance contributions, 
where cos 9 is the angle between the proton and photon in the lab frame, and for the 
CMB spectrum /(e 7 ), the effective threshold is about E p w 10 19 6 eV. The interaction 
length is given by 



and the attenuation length approximately (restricted to two-body final state) by including 
an inelasticity factor K p (E p , e 7 , cos 9) = 1 — E p /E p in the integrand, which is about 1/8 
at the reaction threshold and 1/2 asymptotically. 

The attenuation length sharply decreases with increasing energies up to about 10 20 eV 
due to the rapid increase of the inelasticity factor and cross section. Follows a slower 
decrease up to about 10 20 5 eV due to the inelasticity increase, partly compensated by a 
decline of the cross section. It increases again above 10 21 eV because of the decrease of 
the cross section {K p is saturated). 

The relatively short interaction length ensures that GZK neutrinos are astronomical 
messengers keeping track of the original CR direction. Should there exist a few UHE 
cosmic accelerators located close-by (<Gpc), GZK neutrino detection would allow the 



(17) 



Agzk = (nydpj) 



-l 



3Mpc, 



(18) 
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possibility of pinpointing the location and determining the nature of the most powerful 
cosmic ray accelerators in the universe. 

The GZK neutrino flux has two distinct compo- 
nents, reflected by a double hump structure of the 
expected GZK neutrino flux, illustrated Fig. 4, which 
arises from the kinematic consideration of meson de- 
cay at the higher energies around 10 18 eV (contri- 
bution to v^v^Vp and a bit of v e ) and from the 
neutron decay at the lower energies around 10 16 eV 
(contribution to v e ): while neutrinos from meson de- 
cays carries about 1 /30 of the proton energy (at re- 
action threshold) , neutrinos from neutron decay car- 
ries about 100 times less energy in average than from 
meson decay. 

The exact shape of the GZK neutrino flux de- 
pends on CR composition [53], cosmological UHE 
CR source evolution and the injection spectrum of 
these CR, which is often characterized by a power 
law and the injection cutoff energy. 

These parameters, essential to precisely predict 
the GZK neutrino flux, can be partly constrainedFig. 4. - GZK neutrino flux. Fig. 
from the CR spectral features [54]. However, thefrom [52]. 
energy of CRs is degraded upon propagation and a 
fraction of the CR information is lost and could be 

partly recovered with a characterization of the neutrino flux [55]. We understand the 
necessity of a multi-messenger approach, combining neutrinos and CR information. 

Heavier nuclei are photo-dissociated along their path [56] . Currently, the situation is 
uncertain: the existence of nuclei with E > 10 20 eV is questioned, there may be rising 
evidences for a gradually heavier composition above ~ 10 19 eV [57, 58], compatible with 
an intrinsic limitation of the highest injection energy and/or the GZK mechanism on the 
one hand. On the other hand, a possible correlation of UHE CR sources with the AGN 
distribution by AUGER hints at a light composition [59]. 

Detailed theoretical aspects of this topic are discussed in the lectures by G. Sigl, in 
these lecture proceedings. 

6. Dark matter 

Under certain assumptions, the Minimal Supersymmetric Standard Model framework 
provides a stable weakly interacting dark matter (DM) candidate: the neutralino, a 
self-annihilating thermal relic of the early universe. Its mass, bounded from below by ac- 
celerator constraints and above by theory, lies between 46 GeV [60] up to a few TeV [61]. 
Secondary particles, including v (directly or indirectly), are emitted at a higher rate 
from regions of greater DM density, where gravitationally trapped neutralinos annihilate 
pairwise [62] (neutralinos in the galaxy do pairwise annihilate with a rate which is pro- 
portional the square of the density). The galactic halo or compact objects such as the 
sun seem to be promising regions of such enhanced DM densities and thus for conducting 
dedicated analyses for these signatures, which all share a common approach searching 
for an excess from the directions of these enhanced DM density regions. In the absence 
of any excess over the known atmospheric neutrino background, upper limits on the 
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neutrino-induced muon flux from DM annihilations are obtained. The deeper connection 
to the physics arises with the conversion of the neutrino-induced muon flux limits into 
cross section upper limits [63]: self-annihilation cross sections, velocity averaged (in halo 
analyses), (<jav), and spin-dependent scattering cross sections osd in search for signa- 
tures from the self-annihilating solar neutralinos (assuming equilibrium between capture 
and annihilation rate in the Sun). Assuming neutralinos constitute a sizable fraction of 
the galactic DM density, these analyses have a significant potential to exclude regions 
from the MSSM parameter space which would otherwise remain unconstrained by direct 
search experiments [64] and by 7-ray and CR experiments performing indirect searches 
similar to IceCube [65]. 

6T. Solar WIMPs. - We repeat the old but enlightening derivation of an estimate 
the number of events in a neutrino telescope following [66]: a WIMP halo density p = 
0.4GeVc~ 2 cm~ 3 in corotation with the galaxy v x = 300km/s is necessary to explain 
observed galactic rotation curves, leading to corresponding number density and flux 

500GcV/c 2 \ _ 3 



(19) n x = 8 • 1CT 4 — cm 



, ,'500GeV/c 2 \ 

(20) $ x = nxVx = 2 • 10 4 ( zr -!— ) cm- 2 s" 



Assuming <t x n = (G_fTOn) 2 /to|, we can then derive the solar WIMP capture cross 
section and corresponding rate, 

(21) a sun = ^ = 1.2 . 10 57 x 0.5 • IO- 41 cm 2 

(22) r cap = $ x a sun = 1.2-10 20 s- 1 

for m x = 500GeV/c 2 . 

To finally calculate the rate of solar neutrino of dark matter origin annihilating in the 
center of the sun, we assume a steady state with capture and annihilation rates of WIMPS 
are in equilibrium (the sun has traveled the galaxy multiple times), i.e. r ann = r cap /2. 
The dominant annihilation channel is into weak bosons, each producing muon neutrinos 
with a branching ratio around 10%, XX ~ * WW — > /Uf^, each carrying about half of the 
neutralino energy. The neutrino generation rate is therefore related to the capture rate 
T;/ = Tcap/10 = 1.2 • 10 19 s _1 and the flux at earth is given by 

(23) $„ = ^ =0.5-10- 8 cm- 2 s-\ 
where d = 1 a.u. 

Considering roughly half of the neutrino energy is transferred to the muon (i.e. E„ = 
100 GeV, interaction inelasticity is approximatively 50%), the neutrino-induced muon 
event rate in neutrino telescope with a cross section A = 1 km 2 is therefore given by 

(24) 7V cvcnt = A QvPiceCTv^Rp w 20 yr" 1 , 
where ov^ M = 0.7 • 10~ 36 cm 2 and muon range — 300 m. 
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This is a very small number considering the large atmospheric neutrino background at 
these energies, about 2-3 per year and per square degree and the rather large kinematic 
angle between the neutrino and the induced muon, of order 5 degree. Such a signal could 
eventually be isolated as an excess over background provided a good detector efficiency 
after several years. However, for larger WIMP masses, the potential will slowly degrade 
due to the decreasing statistics partly compensated by the fading atmospheric neutrino 
background. For smaller WIMP masses, it will degrade rapidly due to a number of 
unfavorable factors: increasing atmospheric background, decreasing detection efficiency 
(muon range decrease) and interaction cross section. 

6'2. Neutralino annihilation in the galaxy. - The flux of secondary neutrinos, which 
travel toward the detector from a given direction can be calculated by integrating the 
WIMP annihilation rate along the line of sight 

(25) J(tf) oc J dlp 2 {r{l,$)). 

This approach is model-dependent because the dark matter density profile is not known 
[67, 68, 69]. It is believed to be cuspy at the center and rather flat far from it (where it 
is less model-dependent). Obviously therefore, the high rate would be from the direction 
of the galactic center. Unfortunately, this is a region of high background for the largest 
operating neutrino telescope IceCube located at the South Pole. However the potential is 
not lost, the large background rate being partly balanced by the high expected rate [70]. 
The neutrino flux is [71] 

<»> >>«^<*>^- 

where ((tav) is the thermally averaged annihilation cross section, which is the averaging of 
cross-section times velocity, on the random velocity distribution, the factor 2 because the 
neutralinos self-annihilate and dN v /dE is the neutrino multiplicity on self-annihilation 
of a neutralino pair. The WIMP miracle comes from the fact that the annihilation cross- 
section is of the order of the weak cross-section if we require them to explain the dark 
matter (from freeze out arguments in big bang cosmology). 

7. Neutrino telescope detection methodology 

Neutrino telescopes are instruments optimized for the detection of upward-going neu- 
trino above 10 - 100 GeV, using the earth to filter other particle species (muons), which 
originate in CR showers. The astronomical practicability relies on the quickly fading 
irreducible atmospheric neutrino background with energy compared to expected harder 
ET neutrino fluxes. Current instruments can observe various signatures with limited 
flavor discrimination power and consist of photomultipliers installed in a medium (water 
or ice) transparent to the Cherenkov light emitted by relativistic particles induced by a 
neutrino interaction. The contamination by downward-going atmospheric muons is re- 
duced by placing the detector deeper under ice/water and instrumenting it more densely. 
Typically, the distance between the PMT's will be chosen according to the low energy 
threshold one wants to achieve and the medium optical characteristics. 
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7'1. Neutrino interaction cross sections. - The neutrino-induced muon interaction 
cross section a _, = a(vN — > fJ^X), where Af = |(n +p) is an isoscalar nucleon, above 

a few GeV is [72, 73, 74] 

<-» dMy l:/ ±x ' - ^ (A)' (^)+*^>a-v'» 

written in terms of the Bjorken scaling variables: x = Q 2 /2m-^yE v is the fraction of 
momentum carried by the struck quark (— Q 2 is the 4-momentum transfer between initial 
neutrino and outgoing lepton) and y = 1 — E^/E v is the fraction of neutrino energy 
transferred to the hadronic system (y is called inelasticity factor). Gf is the Fermi 
constant and q(x, Q 2 ) are the quark distribution functions, which have been measured: 

v ^ interaction proceeds with valence quark d v (u„) and sea quark u s (d s ) at leading 
order. 

For illustration, in the range of energy between GeV up to a few TeV (i.e as long as 
E v <C m 2 v /2mN) can be simply expressed as [75] 

(28) dn^/AE^E^E^) « (0.72 + 0.06(£y£„) 2 )10- 38 cm 2 GeV" 1 , 

and with slightly different factor for P M . Therefore (left as a trivial exercice), the total 
cross section 

(29) <T v x{E v )= ^(E^E^dE^ 

Jo d -N 

behaves linearly with energy and the fraction of the energy transferred in average to the 
lepton is 

(30) ^ = 1 - (y) = ±r±- E^{E v ,E»)dE» 

is about 50%. At higher energies, we should account for deviations from the W boson 
propagator. Above 0.1 PeV, the cross section is increasing logarithmically, a vv (E v ) cx 
E® A , and the inelasticity is decreasing (larger energy fraction transferred to the lepton) 
and reach about y = 0.2 above PeV. 

The extrapolation of neutrino cross sections becomes more undcpendable in ultra high 
energy interactions, which emphasize small x < 10~ 4 unprobed regions, together with 
dominant theoretical uncertainty arising from x extrapolation to small values. Within 
the standard model, the current knowledge of the quark structure functions ensures their 
accurate computations up to about 100 PeV [76, 77, 78]. 

Neutrino interaction cross sections with nucleons dominate over the ones with elec- 
trons, due to m p 3> m c , except in the region around E v = m 2 N /2m c w 6.3 PeV, corre- 
sponding to the resonant W boson production 

(31) v e e~ — > W~ — > Dil~ , hadrons 



18 



M. RIBORDY 



(the Glashow resonance): the W boson decay into hadrons (68%) and equally in each 
lepton (total 32%) and the interaction cross section in the resonance region is about 

(32) a{D c e- -> W) w 0.5 • 10~ 30 cm 2 . 

Neutrinos traveling through the earth begins to be absorbed above TeV, due to the rise 
of the cross sections. Modeling the earth density according to the Preliminary Reference 
Earth Model (PREM) [79], the integrated column density X{6) along a chord w.r.t. the 
incidence zenith angle 9 can be calculated. The probability of transmission of neutrinos 
is given by 

(33) Ptr(E v ,6) = c -^N/m P /dip(0 = e -X(6)/A int (E„) 

where the interaction length (in g/cm 2 ) is A int = m p /a u ^. 

The average transmission probability {p iT {E v ))e of an isotropic upward neutrino flux 
(referred to as the shadowing factor S{E V )) is 

(34) S{E V ) = (p tr {E v )) e = J° dcos6e- x ^/ A ^ E ^ 

The case for tau neutrinos is slightly different as they may regenerate when crossing the 
earth [80], see Section 7'2. 

7'2. Event topologies. - All flavors interact with nucleons through neutral (NC) and 
charged current (CC) interactions. Charged current interactions with electrons are for- 
bidden for v^ T . NC interactions and v c CC interaction end up producing shower-like 
signatures referred to as cascades. Muon neutrino CC interactions initiate a shower and 
a muon, resulting into a cascade and long track-like signature. The short lifetime of the 
charged tau lepton and its decay channels will usually only produce a single resolvable 
cascade; only with energies >PeV, the track-like signature may be distinguished (r track 
length of about 50 m/PeV), leaving distinct signatures such as double bang (both the 
initial v T interaction shower and the final shower initiated by the r decay are visible) 
and (inverted) lollipop (only one of the shower is occurring in the detector). 



Muon neutrino. At energies above 10-100 
GeV, neutrino-induced muons is the golden 
channel to search for point sources of ET neu- 
trinos in neutrino telescopes: 

(1) CC interactions result in long muon | 
tracks crossing the detector, collinear with the 
parent neutrino: the average angle between 

and /i from the kinematics is Qv^E^) « 
0.6°/ yjEy/TcY and the long leverage enables 
to achieve resolutions better than a degree as 
shown Fig. 7 and therefore provide a mean for Fi 8- 7. - Estimated IceCube i/„ angular reso- 
astronomy. In contrast, the angular resolution lution - Dash ed: kinematic component, 
of cascade-like events is less favorable. 

(2) Muon neutrinos can be detected when- 
ever the Cherenkov light emitted along the induced muon track is recorded by the sensors. 
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Fig. 5. - Top left: cross sections for CC (solid), NC f M (dashed), CC P M (dashed-dotted) 
and NC £> M (dotted). Top right: inelasticity parameter ycc for muon neutrinos (dashed) and 
anti-neutrinos (dotted). Bottom left: Earth integrated column density X w.r.t. the zenith 
penetration angle. Bottom right: transmission probability of for various zenith angle and 
averaged for an isotropic diffuse flux. 



Combined with the fact that a muon track length surpasses the extension of current de- 
tectors above about 300 GeV, the interaction does not have to necessarily happen within 
the instrumented volume contrary to cascade-like event. This feature significantly boost 
the detection potential of v^. 

For these reasons, we have and will continue to focus mainly on this channel: these 
features combined with the increasing cross section with greater energies maintain some 
constant " dctcctability" for a power law neutrino spectrum &Q/dE u oc E f , with spectral 
index 7 = 2, up to about 10 TeV. Above the " detectability" begins slowly decreasing 
because of significant earth absorption and slower increases (instead of linear) of the 
muon range and cross section. 

The muon energy loss can be parametrized 

d E 

(35) j# = -(a + «U 

where X is the amount of matter traversed by the muon on its way to the detector 
and a w 2.23 MeV/(g/cm 2 ), b w 4.63 • 10~ 6 /(g/cm 2 ) are the nearly energy independent 
standard rock coefficients [81] (of course these coefficients also change with the medium, 
which eventually turns into water or ice) . The critical energy e = a/b w 500 GeV 
marks the transition between continuous and stochastic energy losses. The latter are 
due to interactions of the energetic muon with the (Coulombian field of the) nucleus, 



20 



M. RIBORDY 




Fig. 6. - Shower and track detection modes. Fig. from [83] 



pair producing pairs or disrupting it and bremsstrahlung. Overall, pair production and 
bremsstrahlung processes dominate the muon nuclear interaction. These are of extreme 
catastrophic nature and typically, a 100 TeV muon will loose a significant fraction of its 
energy within only a few stochastic processes. 

Integrating expression (35), we obtain the amount of matter, which is crossed in 
average by a muon of initial energy and final energy E^, 

(36) X(E fi ,E fMl ) = hn^±^. 

R^E^) = X(Q : E^ ) is called the muon range. While at low energy the range 
increases linearly, the range only increases logarithmically when stochastic losses dom- 
inate. However, the range can become very large (the range of a 1 PeV muon in ice is 
i? M w 20 km), much larger than the extension of a detector. 

Showers produced along the muon track through the stochastic processes are produc- 
ing additional light, which makes the precise reconstruction of the direction of the muon 
more difficult and the estimation of the crossing muon energy more precise (continuous 
emission is not providing any hint to the energy of the muon). 

At energies well above PeV, the sky is opening to downward-going muons, due to 
a low atmospheric background, an interesting possibility given earth absorption largely 
prevents the detection of upward-going neutrinos. 



Tau neutrino. Upon v T CC interaction, a large fraction of its energy is transferred to a 
t lepton with a very short lifetime. The r will decay back before loosing a significant 
amount of energy most of the time into a v T neutrino (either directly or through hadronic 
decays) , which keeps in average 20% of the initial energy. Upon NC interaction, E Vt is 
reduced by a factor two in average. Therefore, a v T flux is not absorbed so dramatically 
as the other flavors when penetrating earth at VHE [80] . This feature, called tau neutrino 
regeneration, ensures the experimental possibility of the observation of neutrino point 
sources at VHE, given the equal mixture of all neutrino flavors from cosmic accelerators 
after propagation toward earth (however, the original energy has been degraded and this 
effect only marginally improves the detection potential of VHE neutrinos) . 
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Also, the tau decay branching ratio of into a muon is « 18%. Again, given an 
equal mixture of all neutrino flavors from cosmic accelerators after propagation over long 
baseline distances, the point source search potential at moderate energies, up to 10— 
100 TeV, is slightly enhanced (we are actually only recovering a fraction of the muon 
neutrinos which have been oscillated into v T during propagation). 

7'3. Effective area, event rate and detection potential. - Given a detector with a 
section A and length L, assuming i?„ 3> L, a muon propagating normal to A toward 
the detector will be detected if the neutrino interaction has occurred within a distance 
< i?^ from the instrumented volume. The effective volume can therefore be estimated 
V e s = AR^ (for cascade-like events, in contrast, V c s = AL translating the fact that these 
events are detected if they occur within the instrumented volume). 

The probability of interaction in the effective volume is approximately given by 
p v ^,^(E v ) = R fl ((E t j i0 })/X u (E u ), where \ V {E V ) is the interaction length of a neutrino 
with energy E v (Ai Dt (E u ) = p\ v (E v )) and (E^) is obtained from E v by means of the 
interaction inelasticity discussed above. 

The number of detected events in a time T from a monochromatic neutrino flux 
<& V {E V ) can now be roughly estimated 

(37) N = TAp v ^{E v )* v {E v ). 

The quantity Ap v ^ fJt {E v , (E^)) is called the neutrino effective area A^ s (the neutrino 
flux damping term p tT is neglected in this illustration). 

A less crude approximation of p u ^^ can be made integrating over the differential cross 
section and setting a muon energy detection threshold E^, 

(38) p^(E v ,E,) = -j- f E " dE^(E u ,E^)X(E„E^. 
We obtain for an arbitrary neutrino flux at earth, 

(39) N = TA dE v p v ^{E Ul E^{E v )e- x ^^M. 

Together with the angular resolution function, the effective area is an essential notion 
to describe the potential of a detector. Its knowledge enables to immediately calculate 
the expected number of events from a source, point-like or diffuse. We generically have 

(40) N = T ]T / dE Va A^(E u ,E^e)^(E v ) 

a=C,fl,T V a 

for a point source with flux d^ Ua /AE Ua (for a diffuse source, an integration over the solid 
angle additionally appear). Keep in mind that in order to keep the notation light, we 
have dropped the time dependence of the zenith angle 9 and time integration: in general, 
neutrinos from a point source traveling to a detector will penetrate earth with 0(t), but 
for a detector located at a pole. 
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At this point, however, considering the notion of neutrino effective area only, we 
are missing an essential point (we concentrate on muons) because we can only count 
neutrinos from a given model. But how do we define the search area f2 to compute the 
background in order to optimize the sensitivity when searching for ET neutrinos from a 
neutrino point source and how do we assess whether two models can be distinguished ? A 
neutrino telescope measures the muon energy / incoming direction and not the neutrino 
energy / incoming direction after all. For this, 

(1) The knowledge of the muon effective area A^ S (E^, 9), a geometrically ambiguous 
and ill-defined notion, which can be extracted from 

d$„(E v ) _ r_ Ail . 



(41) J dE v A^(E v ,e)^^ = J dE.A^E^e) 



dE^ 



where the differential muon flux is calculated via the propagation of the muon flux from 
the interaction vertex as 



(42) 



d^^E^O) [,„ m fJ? .d^(E v ) 

— J dE v p det , E v )p tT {E v , 6)p int (E v ) ^ 



PdetiEp, E v ) is the probability density for a muon produced with energy E fl0 (E u ) to 
reach the detector with energy E^, 

(F F ) -dX(E^E, (E u ))/dE, 1 1 



R^ (E V )) \n{l + E^{E u )/e) (E^ + e) 

and set to zero outside of the interval E^ € [0, E^ (E v )]. 

Pint{Ev) = NA<7cc(E v )Rij,(E ft0 ) is the neutrino interaction probability in the vicinity of 
the detector (potentially producing a muon within the reach of the detector) and ctcc 
and a are respectively the charged current and the total muon neutrino cross sections. 
The muon effective area, extracted according to the recipee found in [16] from averaged 
neutrino effective areas found in [82] is shown on the left of Fig. 8. 

(2) The knowledge of the detector energy and angular resolution functions is necessary 
in order to optimize the search area for signal to noise maximization (best detection 
potential) and to calculate the reconstructed differential muon spectrum. 

In [26], we have shown how to extend the analytical treatment above and thus ex- 
tract as well discovery curves, which have been found a posteriori to be in remarkable 
agreement with discovery curves following a full analysis of the IceCube data. They 
are shown on the right of Fig. 8 for various source declinations and bear the following 
meaning: a neutrino flux with spectral index T„ (corresponding to a straight line in the 
log(-E) vs. E 2 log( Flux) representation) tangent to one of the spectral discovery curve is 
at the limit of discovery. This representation of the IceCube discovery limit (s) is useful, 
because it enables 

1. The estimate of the sensitivity for arbitrary assumptions about the slope of the 
neutrino spectrum (i.e. the normalization of the minimal detectable spectrum at a 
given reference energy) 

2. The estimate of the neutrino energy range contributing most significantly to the 
source signal (the energy at which the minimal detectable source spectrum grazes 
the discovery curve). 
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Fig. 8. - Left: Full lines are the extracted muon effective areas for IceCube, when in construction 
(22 strings) and completed (80 strings). The dashed curve is the muon effective area derived 
in [51]. Fig. from [16]. Right: IceCube point source reach after three years of exposure. The 
tangent to one of these curves intersecting the y-axis at 1 GeV defines the muon neutrino flux 
normalization for the corresponding declination. Fig. from [26]. 



The dependence of the discovery curves on the source location in the sky is reduced 
to the source declination 5, at a polar location. We notice a weak dependence of the 
discovery limit on the source declination in the case of soft neutrino sources (T„ ^> 2): 
Most events contributing to the source signal have relatively low energies, at which earth 
is transparent for neutrinos. For sources with hard spectra, r„ <C 2, the discovery 
potential is strongly affected by the increase of the source declination. The rise of the 
earth absorption probability with energy obscures the source signal, contributed mostly 
by high energy neutrinos. The effect is most dramatic for the hardest neutrino spectra, 
r„ ~ 1. In this case, the minimal detectable flux for the sources at the declination 
S = 75° is about 1.5 orders of magnitude higher than from sources at 5 = 15°. 

Experimentally, all of the complications above vanish and the potential of the detector 
and its power to reject models or delineate concurrent model comes straightforward from 
the analysis, relying partly on the detector Monte Carlo (mainly for the simulation of 
the signal). 

7'4. Event reconstruction. - Tracks or showers are reconstructed using a probability 
density distribution (PDF), which depends on measured time and location of the detected 
Cherenkov photons (hits). 

Consider a Cherenkov light source (e.g. a track or a shower), which can be parametrized 
by giving a vertex q, incidence angles 9, 4> and its energy E, that is rid.o.f = 6. For a 
source with a different topology, e.g. an extended source, the parameters should be 
defined accordingly and rid.o.f acquires a corresponding value. 

The reconstruction of the parameters consists in maximizing the log likelihood ex- 
pression for an arbitrary hit time series for each receiver [83], 

(44) ta £_^^(f;( hlA(+1 . Aj ) + Y. ^i) 

d.o.f = 1 {j\N j= 0} J 

where j S {1, .., M } is the hit receiver index, nhit = Y] j Nj is total number of hits in the 
event, p* N . = p*(tx, ...t^) oc Ili^iPi ^ s the PDF of detecting Nj photons at receiver j 
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and at times {t\, ...t^} (pj is the PDF of detecting a photon at receiver j and at time 
ti) and fi is the probability to detect n photons at receiver j, assumed to obey Poisson 
law — e~ M /i"/n!, where fi = /jP is the mean number of detected photons, pj's and 
yu's depend implicitly on the relative orientation and distance between the source and 
the receiver, the optical properties of the medium, the receiver efficiency, etc. The yu's 
additionally depend on the source intensity (related to the source energy). Some detail 
on the specific expressions of these PDF's as used in IceCube can be found in [84]. The 
last term represents the information from receivers which have not detected any photons. 

We notice that eq. (44) can be split in two parts depending separately on energy E 
and directional {q, 8, (f)} parameters. This leads to introduce the reduced log likelihood 
formulas: 



1 M 

(45) In £ dir = V ]np> N 

"hit - n d . . f ^ 

and 

(46) ln£ g= 1 f (f>4+ E 

where n<j. .f take the corresponding values (respectively 5 and 1). 

Expressions (45) and (46) can be used separately to reconstruct the corresponding 
subset of parameters in a first approximation, when calculation speed or convergence is a 
concern. This will consequently be accompanied by a degradation in the precision of the 
reconstructed parameters and best results are obtained by using expression (44), which 
more fully exploits the available information. 

A more tractable, hybrid, so called MPE (multi photo-electrons), reconstruction, 
which leads to very good reconstruction results is based on ln£dir, but with a modified 
PDF fP N . . This PDF does not use explicitly all hit times but only the first one t\ and 
the total number of hits Nj. It can be simply derived from the above and one obtains 



(47) 



where P[ is the cumulative of p\ . 

Performance of the reconstruction in IceCube will be explicitly illustrated in a subse- 
quent Section 10'3. 

Generalization. The log likelihood expression can be formulated generically for the case 
of composite events [85], i.e. events whose mixed detector response is due to the jux- 
taposition of more than one Cherenkov light source. Such a reconstruction is useful to 
reconstruct an arbitrary event topology and to favor or discard a given event topology 
hypothesis (by contrasting In C of the various hypotheses) . The likelihood minimization 
becomes slightly more difficult due to the increase of parameters and the fact that ln£ 
cannot be split into the two directional and energy components, which are entangled. 
Composite topologies will result for instance from 
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1. Down-going uncorrelated muons (originating from distinct atmospheric showers) 
track events. We expect a rate of a few per second at 2 km depth in cubic km-scale 
detectors; 

2. Events from muons undergoing catastrophic energy losses leave a signature resulting 
from the superposition of showers and a track; 

3. Neutrinos converting in the detector effective volume (a shower and a track, possibly 
finite); 

4. Exotic channels (stau pairs, micro black hole evaporation) may result in nearly 
parallel track events [86, 87]; 

5. PcV tau double bang and lollipop events. 

7'5. Systematic uncertainties. - Besides intrinsic difficulties (e.g. kinematics of the in- 
teraction) and statistical limitations (ET flux versus atmospheric neutrino background) , 
various source of systematics uncertainties contribute to the adversity faced by anal- 
yses, which result from our limited knowledge of the instrumental and physical input 
parameters: 

- Bioluminescence in water, optical properties in the ice: dependence of absorption 
and scattering profile and length w.r.t. the depth and refrozen ice holes; 

- Optical module absolute sensitivities and their exact locations and orientations; 

- Interaction cross sections (rising uncertainties at UHE); muon propagation; 

- Atmospheric background spectral shape and absolute normalization (cosmic ray 
spectra, composition and subsequent muon and neutrino production). 

These uncertainties cannot be reduced by increasing statistics, at least not in a straight- 
forward manner: some are in principle reducible indirectly with increasing statistics, 
such as the ice properties, the difficulty residing in the decoupling of entangled uncer- 
tainties. 

8. Neutrino telescopes: past and present 

Greisen posits the concept of neutrino telescopes in 1960 [88], "...we propose a large 
Cherenkov counter, about 15 m in diameter, ..." "...about 3000 tons of inexpensive liquid. 
..." "... from the Crab nebula the neutrino energy emission is expected to be three times 
the rate of energy dissipation by the electrons, ..." and concludes "Fanciful though this 
proposal seems, we suspect that within the next decade, cosmic ray neutrino detection will 
become one of the tools of both physics and astronomy. " 

Reines discusses simultaneously the arguments, which remain valid nowadays [89]: 
"... [cosmic neutrinos] propagate essentially unchanged in direction and energy ... and 
so carry information which may be unique in character. For example, cosmic neutrinos 
can reach us from other galaxies whereas the charged cosmic ray primaries reaching us 
may be largely constrained by the galactic magnetic field ..." and notice, due to the lack 
of understanding of the origin and propagation of the charged cosmic rays: "... the 
cosmic neutrino flux can not be usefully predicted. " and add "The situation is somewhat 
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simpler in the case of cosmic-ray neutrinos{ A ) : they are both more predictable and of less 
intrinsic interest. " We appreciate now the irony of this statement. 

8'1. First generation. - In 1976, the precursor DUMAND (Deep Underwater Muon 
And Neutrino Detector) project was initiated [90]. The early sketch had the ambition 
of a cubic kilometer telescope, but was downscaled several times. A prototype string 
(DUMAND-I) was successfully operated from a ship off Hawaii in 1987 at a 4 km depth. 
In 1993, elements for a 9-string detector array DUMAND-II consisting of 216 photomul- 
tiplier tubes and corresponding front end electronics encased in pressure spheres were 
essentially ready. Deployment of the infrastructure and a string began in December 1993. 
Unfortunately, the leaking of an electrical penetrator and the subsequent failure of the 
string controller electronics lead to the interruption of the observations after 10 hours of 
operation. The project was eventually canceled in 1994 due to a lack of funding. 

Parallel to the DUMAND project, the Baikal neutrino telescope effort started in the 
early 80's [91] and was the first to observe high energy neutrinos [92] in a four string 
configuration in 1996. The detector located 1.1 km underwater at the Baikal Lake in 
Russian was expanded into the NT200 configuration in 1998 and then into NT200+ 
configuration in 2007 with the addition of three sparsely instrumented outer strings [93] . 
Maintenance and deployment take place from the frozen surface of the lake in winter. 

The realization of the AMANDA (Antarctic Muon And Neutrino Detector Array) 
project was successful after overcoming a certain number of difficulties. In contrast 
to DUMAND and Baikal, AMANDA was based on an alternative concept, using the 
Antarctic ice cap rather than water as the deep underground detector medium and tar- 
get, escaping some of the limiting water backgrounds (bioluminescence and seawater 40 K 
decay). A first string equipped with four optical modules (pressure sphere identical to 
DUMAND) was deployed in 1992 and successfully operated at 800 m depth. In 1994, four 
additional strings with 20 optical modules each were installed (AMANDA- A) , including 
optical fibers for laser calibration. It quickly became clear that the characteristics of 
the ice at this depth were not appropriate because of the contamination with air bubble 
scattering centers: the Cherenkov light effective scattering length lying at values about 
20 cm, i.e. light was heavily delayed when propagating, prevented any reliable recon- 
struction of the incidence muon directions and therefore efficient searches for neutrinos 
to be conducted [94]. This experience nevertheless confirmed the viability of the concept. 

Ice cores showed that bubble size and density were decreasing drastically deeper in 
the ice. The prospects of a deep detector would be further improved considering the 
additional advantage of a reduction of the atmospheric muon background. In 1997, the 
AMANDA-B10 detector was deployed and operating [95, 96] and it observed atmospheric 
neutrinos [97]. In 2000, an upgrade phase started and nine more strings were deployed 
(AMANDA-II). Strings are located at depths between 1500 and 2000 m. Most notably, 
AMANDA-II has produced very stringent limits on diffuse neutrino fluxes [98, 99] and 
neutrino flux from point sources [100] and accurately measured the atmospheric neutrino 
spectrum [101, 102]. The absence of evidence for ET neutrinos has provided constraints 
for the future of high energy neutrino astronomy, setting detector scales to cubic kilometer 
instrumented volumes. 

A recent review [103] presents the initial pioneering neutrino astronomy phase in 
greater detail. 



( 4 ) note: atmospheric neutrinos 
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Fig. 9. - The AMANDA neutrino telescope. 



8'2. Present generation. - There are two principal neutrino telescopes in operation: 
Antares in the Mediterranean and the IceCube neutrino observatory, which signals the 
emergence of a new class of gigantic detectors dedicated to the observation of the high 
energy (HE) neutrino sky. The vast research program is found in [104, 105, 106, 107, 
108, 109, 110]. 

IceCube. The IceCube observatory was completed in December 2010, after a very suc- 
cessful deployment over the past five years [111, 112]. The in-ice array is equipped with 
86 strings, a nominal string spacing of 125 m Each string is equipped with 60 optical 
modules, enclosing a photomultiplier and digitization / communication / time stamping 
boards [113, 114]. The design, construction and deployment of IceCube has largely prof- 
ited of the gained expertise from the AMANDA enterprise at the South Pole [96, 97], 
which operated between 1997 and 2009 and provided a proof of principle [115]. 

In the center of the in-ice array, 8 strings more densely instrumented form together 
with the neighboring strings a dense 20 Megaton inner core, enhancing the IceCube 
detection capabilities toward lower energies (the data are enriched with lower energy 
neutrinos) and potentially enabling IceCube to explore the Southern neutrino sky. The 
design is meeting performance expectations [116] and IceCube is taking data at a rate of 
>2 kHz, dominated by CR muons. 

On the ground surface, on top of in-ice IceCube array, the IceTop air shower array 
instruments a square km. IceTop is dedicated to composition studies of the CR spectrum 
around the knee and above [117]. IceTop enables the detection of air showers at energies 
above ~1 PeV. As a part of IceCube, it also plays an important role of facilitating 
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calibration of events triggered by the in-ice array [118] and background rejection of very 
high energy events. 

At the deployment depths, the secondary muon intensity originating from cosmic air 
showers in the atmosphere is reduced by a large factor. This flux remains the major 
experimental background for these detectors. The soft atmospheric neutrino background 
constitutes an irreducible background for ET neutrino searches, which can only be com- 
pensated by expected harder signal energy spectra. 

IceCube's main goals are the unambiguous identification of the first galactic and extra- 
galactic cosmic ray (CR) accelerators with the detection of HE neutrinos from point 
sources and the divulgation of the nature of dark matter (DM) through the observation 
of a secondary neutrino flux from annihilating DM in our galaxy. 

IceCube 



ceCiibe Lab 




Fig. 10. - Left: IceCube. Right: comparison of the IceCube point source flux discovery potential 
for one year with IceCube in its 40 string configuration and after ten years with the completed 
array for different source declinations. 

The rate of increase of the IceCube discovery potential has been greatest during the 
deployment phase and will continue to be high until about 3 years of data have been 
analyzed. A noticeable that a jump by an order of magnitude of the discovery potential 
will occur within the next 10 years (for a 20-fold increase in the statistics), see Fig. 10, 
right, compared to current results [119]. In the case of point source discovery during this 
period, an era of precision for HE neutrino astroparticle physics will be initiated and 
the prospects will include the characterization of the detected neutrino fluxes (i.e. the 
energy spectra and their time-dependence). 

On site, prospecting activities for the detection of UHE neutrinos are taking place 
as well, relying on the alternate radio Cherenkov and acoustic signatures accompanying 
UHE neutrino interactions, see Section 11. 

Antares. Antares operates in the Mediterranean See since May 2008 at 2.5 km depth and 
consists of 900 PMTs distributed on 25 storeys and 12 lines (3 PMT per storey), placed 
60-80 m apart, thus instrumenting a volume of 0.025 km 3 . Not located at a Pole, the 
field of view is continually changing and exposure is not uniform but complementary to 
IceCube. It consists of 13 strings half a kilometer long, deployed at a depth of 2.5 km. 
An excellent muon angular resolution of 0.1° can be achieved due to the low scattering 
of Cherenkov light. In this non sterile environment, one has nevertheless to cope with 
currents, organic lightning, dirt deposit on the optical modules and with potassium decay, 
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leading to 100 kHz PMT noise rate, necessitating stages to be equipped with 3 modules 
operated in coincidence. 

9. Supernova neutrino detection methodology 

The core collapse supernova energy release is mainly directed into the neutrino com- 
ponent. Estimates posit the frequency of these events in the Milky Way of about 2 ± 1 
per century, so that expectations are optimistically 25% probability of such an event in 
the next decade [120]. Such an event would be of great significance, not only constrain- 
ing the dynamics of the core collapse itself but potentially granting access to neutrino 
hierarchy (remember that the sign of difference of mass eigenstates is accessible through 
matter oscillations and scales with #13). 

The main phases of 10-100 MeV neutrino emission are successively 

1. O(10 ms) deleptonization phase, a burst of electron neutrinos; 

2. Few 100's ms accretion phase; 

3. Cooling phase, with roughly all flavor energy equipartition. 

Far more detail on the physics of supernova can be found in the dedicated lectures by 
G. Raffelt in these school proceedings. 

9'1. Neutrino telescope as supernova detector. - Neutrino telescopes operating in the 
quiet ice medium can be utilized as instruments for supernova detection, due to the very 
large neutrino flux from galactic supernova at earth and moderate optical module noise, 
of order r — 0.58 kHz, this even while the cross section are very small and low photon 
yield signature of individual event is inappropriate for these types of instruments. 

The traditional detection principle is based on the sudden increase of the detector 
hit noise rate, as originally proposed [121, 122, 123] and IceCube can monitor the whole 
Milky Way [124]. The main detection channel is the interaction of anti-electron neutrinos 
with hydrogen protons, 9 c p — > e + n. The contribution of other channels is suppressed 
due to the large oxygen 16 binding energy, the rare presence of oxygen and hydrogen 
isotopes and the small electronic cross sections. A significant fraction of hits (about 10%) 
arise from neutral current cross section with oxygen. 

Mediterranean telescopes cannot be effectively used the same way due to the large 40 K 
radioactivity and microscopic bioluminescence (both components accounting for about 
100 kHz optical module noise rate in quiet conditions). 

The brief and powerful electron neutrino neutronization emission is of great impor- 
tance given its connection to astrophysics and particle physics and can be studied via 
electron recoil in elastic collisions and to some extent through oxygen charged current 
reaction. The burst luminosity also moderately depends of the progenitor mass, granting 
potentially access to the distance of an obscured core collapse supernova located at 10 
kpc with 5% precision [125]. 

Let's consider a benchmark supernova at d = 10 kpc releasing E$n = 2 x 10 59 MeV 
in neutrinos within about T = 3 seconds (a nearly optimal integration time in a simple 
analysis to assess the detection potential in SF model [126]), we can calculate the time- 
averaged D e flux, 



(48) 
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where factor 1/6 reflects flavor equipartition. 

In the energy range of interest, the interaction cross section is quadratically growing 
with energy, about <7p oP s=s 10~ 41 cm 2 at 12.5 MeV. We obtain the interaction rate in a 
gigaton detector volume V with density p, 

(49) pV^ e a PeP N a /9 w 0.7 • 10 8 s" 1 

Note that it would be more correct to use in this case a cross section weighted on the 
differential neutrino flux, 

(50) J &E v j AE c d ^(E v ,E c ) ^{E v , d) 

in place $p c crp oP , but we drop this for simplicity in this back of the envelope calculation. 
It would also turn out that the average positron energy is significantly higher than the 
average neutrino energy as shown in Fig. 11. 

Positrons from 12.5 MeV neutrinos will roughly have energy E c = E v — (Am pn + m e ), 
where Am pn + m e ss 1.8 MeV, i.e. emit about N(E C ) = 2000 Cherenkov photons between 
300 and 650 nm along a 5 cm track. Due to multiple scattering, the Cherenkov photon 
emission will be quite uniform over tt sr. 

The optical module in IceCube have a diameter of about 2R = 30 cm and a photo- 
detection efficiency of about £pmt = 7% (photo-multiplier tube quantum efficiency and 
glass), once integrated between 300 and 650 nm for a Cherenkov spectrum. It is simple 
to estimate the effective volume of the IceCube detector composed of about Nq M = 5000 
modules for the detection of uniformly interacting neutrinos in the detector volume, 

r 2 dr J dcos0/(>fc, r, 0,E C ), 

where /(> k, r, 0) is the probability density of recording at least k hits from a positron at 
distance r from an optical module and at an angle (optical modules are not uniformly 
sensitive, so we define 9 as the angle between the optical module axis and the positron 
with the optical module at the origin of the coordinate system). 
For k = 1, we have 

(52) /(> 1, r, 6, E c ) « /(l, r, 6, E c ) « e PM TN(E c ) e(r, 9) e~ r ^ 

where e(r,6) is the fraction of the 47rsr solid angle occupied by the sensitive surface of 
the module, viewed under angle 9 and at distance r and A a b s is the medium absorption 
length of the Cherenkov photons. is slightly above 2 megatons, i.e. about 0.2% 
of the instrumented volume, or a hit rate of 200 kHz, well above the nominal standard 
deviation of the detector noise rate ^/V M( r /Hz) Hz w 1.7 kHz. Integrated over the 
optimal time window T, we may therefore expect a 5 standard deviation excess up to 
distances of about 65 kpc. Detailed analysis presented in [124] confirms this picture. 

The traditional method is therefore very effective for the detection of galactic super- 
nova and IceCube complements other operating neutrino detectors, e.g. Super-K (few 
thousands detected neutrinos for this benchmark supernova), useful when an instrument 
is in maintenance, to cross calibrate, triangulate the SN incidence direction (IceCube 
resolution is 2 ms) or to study earth matter effect. 
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Fig. 11. - From left to right. SF model neutrino emission spectrum (thin, dotted) and positron 
spectra: produced (thin, dot-dashed), producing single hit (solid), double hit (dashed), 1 + 1 
(dot-dashed) and 2 + 1 (dotted) signatures for IC. 



Assuming an emission spectrum, this method enables the extraction of information 
related to the supernova luminosity as well, once the supernova distance is known. How- 
ever, the neutrino flux and average energy cannot be disentangled. Moreover, in the 
case the supernova escaped optical detection, the luminosity will remain unknown. It is 
therefore interesting to consider alternative detection channels. 

9'2. SN detection based on coincident hit rate. — Coincident hits are hits at a single 
or neighboring optical modules, occurring within a very short time window, of about 
< 100 ns and originating from a single positron. The ratio of coincident to single hit 
rates increases with the average positron energy: while the single hit probability is rising 
linearly with the positron energy (i.e. the positron track length and number of emitted 
Cherenkov photons), the probability for two correlated hits is rising with the square of 
the positron energy. 

Therefore this ratio constitutes a new observable for the measurement of the average 
energy of the supernova neutrinos. Combining it with the single hit rate, it enables to 
disentangle neutrino flux and average energy, in contrast to the traditional supernova 
detection method. 

Let's discuss here only the case of two coincident hits at a single module for sim- 
plified treatment and Sf neutrino emission spectral shape model: the average energy 
of the positrons, about 17MeV, is significantly larger than the neutrino energy, about 
w 12.5 MeV, due to the rising neutrino interaction cross section with energy. Average 
energies of the positrons resulting in a single and double hits are respectively about 20 
and 22 MeV, illustrated in Fig. 11. Using these numbers and eq. (51), the ratio of single 
to double hits can be calculated and one obtains about 1%. This ratio sensitively depends 
on the average energy and few percent accuracy would be obtained for the benchmark 
supernova. 

9'3. Backgrounds and other difficulties in SN analyses. - The analysis of supernova 
neutrinos has to cope with the atmospheric muon background, which contribute about 
15 Hz to each optical module hit rate in average. This background is also at the origin of 
correlated hits in neighboring optical modules and can be reduced by removing external 
detector layers (atmospheric muons dying in the superficial detector layers and producing 
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hits insufficient to fulfill the event trigger condition). Seasonal variations of module hit 
rates, due to seasonal change of the temperature and pressure in the atmosphere, also 
affect the analysis and outcome in term of significance of an observation and must be 
carefully accounted for. Sensor correlated hits at second and milli-second time scales, due 
respectively to photomultiplicr aftcrpulse and radioactive decay chains may eventually 
also reduce the coincident hit method potential. 

Calculation details, rejection of correlated hits from muon background, results from 
simulation and prospects with future IceCube extensions are found in [127]. 

10. Point source search methodology 

lO'l. The multi-messenger approach. - The existence of a correlation between the 
energy released by an astronomical source in the form of electromagnetic (EM) radiation, 
CRs and neutrino flux is at the origin of the multi-messenger (MM) paradigm, supported 
by convincing phcnomcnological arguments and numerous models of source activity. 

The most basic expression of analyses performed in the MM context is the search for 
a significant excess of neutrinos (over the atmospheric neutrino background) from the 
direction of point source candidates selected on the basis of their MWL spectra along 
with the appropriate phenomenological arguments. This is diametrically opposed to a 
random search for neutrino point sources (e.g. a whole sky search), which has no direct 
connection to phenomenology (with IceCube, it has been demonstrated that the analysis 
method is robust with a change of the assumed spectral index as we will see later). The 
MM approach augments the significance of observation, however, this is true only if the 
actual neutrino sources are among the selected targets: the approach relies on the validity 
of phcnomcnological arguments; moreover, the existence of neutrino sources not emitting 
in the electromagnetic spectrum cannot be excluded. 

10'2. Point source search strategies. - Point source searches are analyses looking for a 
statistical excess originating in narrow sky regions. Best current sensitivities are rather 
uniform w.r.t. the declination and notably below the intrinsic TeV radiation strength 
of e.g. Mrk 501 during the 1997 flaring state (assuming "f/v ss 1), These searches 
are optimized by taking advantage of the experimentally observed off-source detector 
response, which defines the atmospheric neutrino background i.e., at first order, sources 
of systematic errors can be neglected. 

Steady point source searches, on the whole sky or based on catalogs, are conducted 
within the HE neutrino telescope community [119, 128, 129, 100, 130]. Concurrent 
searches, based on MWL spectra, introduce an optimization based on the spectral char- 
acteristics and/or the time dependence of the photon emission in a certain energy band 
reflected in the neutrino emssion, for instance for GRB's, periodic galactic systems, flar- 
ing objects such as AGN's [131, 132, 133, 134, 135, 136]. 

10'3. Whole sky steady point source search. - We discuss now the generic search 
methodology for localized excess in the sky in greater detail. We start with the very 
basic approach, so called "binned" search and restrict to the very specific South Pole 
location for simplicity, as the background density is a function of declination 5 only: the 
zenith angle of the incoming neutrino direction is related to <5 independent of time 
and the right ascension a drops assuming a cylindrical detector symmetry and steady 
emission. 

Consider a circular or square bin of area a around a source located at (6, a). The 
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bin extension is chosen in relation to the point spread function Q Vl i(E) for the entire 
neutrino sample at this declination. Given N neutrino candidate events in the sample in 
the S declination band of area A, the bin background expectation is 

(53) Mb = Np = N(a/A). 

The significance of an observation of n Q b s events in the search bin in the absence of 
signal is obtained using the cumulative distribution function (CDF) Ib of the binomial 
probability Ps(n bs; N, /ib)- In point source searches, the significance of an observation 
can be set without relying on the Monte Carlo. Obviously this is not the case when it 
comes in setting a flux or flux limit: a detector response model is necessary. In practice, 
the event candidate selection is made with the help of the Monte Carlo, in order to adjust 
the analysis to a specific assumption on the source spectrum (usually taken as a power 
law spectrum with 7 = 2 extending to very high energies, eventually dominating the 
background of atmospheric neutrinos) 

In order to avoid any bias in the selection of the neutrino candidates, the events 
selection optimization is performed "off-source", using randomized right ascension a, 
which is, still restricting the approach to the South Pole location, equivalent to scrambling 
the recorded event timestamps. 

Concept of upper limit and sensitivity and discovery flux. Once one disposes of an actual 
measurement, assumed to be below the discovery threshold arbitrarily set to something 
like a 5a significance of non zero observation, the upper limit at some confidence level 
(C.L.) can be set (usually something like 90%). The meaning of the 90% C.L. upper 
limit is the following: for a given experimental outcome n b s , the true unknown signal 
lies below the 90% flux upper limit $90 (or p, 90 in terms of number of events) in 90% of 
the experiments [137]. The construction of upper limits presented in the literature often 
relies on the Neyman [138] or Feldman-Cousin [139, 140] prescriptions. The conversion of 
/igo into $90 necessitates the Monte Carlo as already stated and must assume a spectrum. 

Before the actual measurement takes place (i.e. the experimental data analysis), con- 
current analyses must be tuned and performance compared. A useful concept is the 90% 
average upper limit or sensitivity. As the name indicates, the sensitivity is the upper 
limit averaged for all possible outcomes n b s of an ensemble of experiments in the absence 
of true signal, 

00 n obs 

(54) M9o(Mb) = Ai9o(«obs, Mb) -^-j-c^ b 

n ^obs- 

The power for discovery of an experiment after a certain operating time is another 
useful analysis outcome. The literature sometimes mentions the 5cr discovery flux in 
50% of the experiments. In the simple case of a binned search, this can be built 
easily: the binomial PDF defines the minimum n bs,5t7 for a detection at the 5cr C.L. 
given fi h . The discovery flux is set by finding the right amount of signal /li s to add to the 
binomial PDF 

(55) fs(Mb) -> fs(Mb + Ms) 

necessary to ensure an experimental outcome n Q b s > rc bs,5<T in 50% of the experiments. 
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This binned search is appropriate to obtain the results analytically, i.e. building the 
test statistics (TS) Ib and Ib and (x s can be extracted analytically 

(56) [M s | {Js(/Xb + A* S ,"obs,5 CT ) = -j}, 
where 

(57) n bs,5<7 = "obs | {hi^b, n ohs ) = ^erfc(5/V2)} (@ 5a C.L.). 

In practice, and in order to introduce the concept once we will consider the more 
powerful "unbinncd" search methodology (which, as a drawback, relies heavily on the 
Monte Carlo), the construction of the TS illustrated Fig. 12 proceeds according to the 
following steps: 



7 B(»W W, Mbg) 




Fig. 12. - Illustration of the building of the test statistics in the simplified case of a binned 
search. 



1. By means of experimental and Monte Carlo data, the TS is generated by scrambling 
the sky maps of experimental data million of times including some variable amount 
of Monte Carlo signal according to the assumptions of source location and spectrum; 

2. For each map, the TS distribution is filled with n Q b s in the search bin around (5, a). 

The discovery flux is found by comparing the null signal CDF with the set of CDF's 
with non vanishing signal. The discovery flux leads to rt ODS > "obs,5cr in 50% of the 
experiments, where n b s ,5o- is the number of events in the absence of signal that would 
be attained by a fraction 1/1.7 • 10 6 of the experiments. 

Generalization: "unbinned" search. We discuss now the straightforward generalization 
of the methodology to "unbinned" point source searches [141]. The additional power 
originates in the non binary underlying methodology: in the binned method, no maximum 
profit was taken of the shapes of the angular PSF and energy PDF, i.e. we only minimally 



METHODS AND PROBLEMS IN NEUTRINO OBSERVATORIES 



35 



exploited our knowledge of the detector and physics (kinematics, muon energy at the 
detector). We now apply a maximum likelihood fit to the data to determine the relative 
signal contribution. 

We assume a data sample, which consists of N neutrino candidates, as a mixture 
of two components, signal and background. The probability density of the ith event 
measured at location ((5j,aj) with energy Ei, assuming a signal source with spectral 
index 7 located at (<5, a) is 



(58) 



^ + ^ - 



where the signal and background PDF's are 



(59) 
(60) 



Si = M s {{8i, a t } I {5, a}, a,) x £(Ei | 7, Si) 
Bi = 7V at m(^) x £ (Ei I atm, S t ) 



Here Af s is the probability that an event originates from the source and modeled as a 2-D 
gaussian of width er^. <jj is the angular resolution obtained for each event individually 
from the error on the parameters from the maximum likelihood reconstruction fit( 5 ); J\f a 
is the angular PDF for the atmospheric background and is flat w.r.t. a at the South 
Pole; £ describes the probability of the reconstructed energy Ei . 

Fig. 13, left, shows the PDF of the muon energy estimator for various simulated signal 
power laws and for experimental data, while Fig. 13, right, shows for simulated data the 
correlation between <Ji and the actual track reconstruction error. 
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Fig. 13. - Estimators. Figs from [119]. 



The parameters n s and 7 are obtained from a maximum likelihood on all events 
i € {1, N} using the PDF from eq. (58), 



(61) 



JV 

£(n.,7) = Il[ 



—St 
N 1 



(1 - ~)^, 



( 5 ) More precisely, at is obtained from a faster algorithm (not from the error processors available 
in the MINUIT package), the result from a paraboloid fit of the likelihood landscape space of the 
reconstructed incoming direction of the event (0, <f>) is an estimator of the angular uncertainty. 
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The TS distribution, built with the millions of scrambled sky maps, is commonly defined 
against the null hypothesis n s — 0, 



(62) 



A = -2 In 



C{n s = 0) 



In order to assess the significance of an observation, the following steps must be taken: 

1. Build the TS distribution. Compared to the binned search, it simply consists into 
the substitution of n bs with the two-step process of maximization and calculation 
of A; 

2. The final number, the p-value associated to A is the fraction of scrambled data sets 
leading to higher TS values than A. A defines the discovery flux, and used in the 
Feldman-Cousin construct for instance, the sensitivity and upper limit. 

These are illustrated on Fig. 14, left: we see the TS distributions of pure background 
and background + various signals; at A s» 24 — 25, the residual fraction of experiments 
leading to larger A is at the 5cr level. The signal TS shows that a source with a strength 
corresponding to a number of selected events between 8 and 12 is necessary to have 50% 
of the experiments with A > 24 — 25. 

Fig. 14, right, demonstrates that this formalism enables the study of the resolution 
on 7 and the number of selected source events required for discovery for various spectra, 
7 = 1.5, 2, 3. 
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Fig. 14. - Test statistics in the case of an unbinned search and potential of the method (see 
text). Figs from [119]. 



A further complication arises, related to the statistical trial factor. If a search looks 
for evidence of an excess from a small number of predefined and well separated sources, 
the trial factor is the number of sources in this catalog. When searching in the whole 
sky for a point source excess, it is not so well defined: for a given monochromatic source, 
the "effective" number of sources would roughly correspond to An/a(E), where a is the 
subtended area from the PSF width. In practice, the unbinned analysis is repeated on 
the global sky, filling the TS distribution with the sky location corresponding to the best 
p-value. One obtains the post-trial p-value. 

Fig. 15 shows the IC-40 sky in terms of 90% upper limits. One notices that they are 
best for sources at positive declination and are degraded by 10-100 in the Southern sky 
for obvious reasons. 

The unbinned search methodology can be extended in a very straightforward fashion, 
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1. For transient and periodic sources (for instance GRB, AGN burst, binary systems), 
by simply adding the corresponding time-dependent factor in the signal PDF. 

2. To stacked searches, which consist in cumulating the signal from similar sources 
with different weights. E.g. the weight may depend on the source distance or on 
the source acceptance (different background at different declination at the South 
Pole, different exposure in ANTARES). 

3. To extended point source searches, by replacing the PSF with the PSF convoluted 
with the source distribution: PSF — » PSF ® source distribution. 

11. UHE neutrino search with alternative detection techniques 

Typical predictions for the cosmogenic neutrino flux, i.e. neutrinos from photo-pion 
production, is of the order EdN/dE ~ 10~ 17 s _1 cm _2 sr -1 at E = 10 18 eV. In ice, this 
flux translates into one neutrino interaction / km 3 every few years. As we mentioned 
previously, however, the prediction strongly depends on the primary cosmic ray com- 
position among other parameters and the rate could be lower by one or two orders of 
magnitude. 

The characterization of the GZK neutrino flux spectrum, and thus the recovery of the 
degraded information carried by UHE CR's, would permit the delineation of cosmological 
source evolution scenarios from source injection spectrum characteristics. To fulfill this 
goal, however, the GZK event detection rate should be vastly increased. Therefore, 
a much larger volume should be instrumented with an appropriate technology for the 
detection of ultra high energy neutrino interactions. 

As the applicability of the optical detection technique reaches its limits due to < 100 m 
attenuation length, two alternative detection methods were proposed for operation on 
large arrays, following signatures first discussed by Askaryan [142, 143, 144]: an inter- 
acting neutrino emits a coherent Cherenkov pulse in the range of 0.1-1 GHz [145] close 
to the shower axis and a thin thermoacoustic pancake normal to the shower axis [146]. 

The bremsstrahlung and pair production cross sections starts decreasing quickly above 
O(PeV) strongly affecting the longitudinal profile of the shower above EeV energies, 
which then consists of disparate subshowers initiated by these UHE daughters. This 
feature, known as the LPM effect (Landau-Pomeranchuk-Migdal) [167, 168] and funda- 
mentally related to the non point-like nature of the interactions, reduces in particular 
the prospects of the acoustic detection technique (higher energy threshold). 

The radio detection technique exploited by several past, operating or projected in- 
struments observe the direct [147, 148, 149, 150] or the refracted radio Cherenkov pulse 
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either from the air [151] or from space [152] following an interaction near the earth's 
surface, or they observe the signal following an interaction in the lunar regolith [153]. 

Acoustic sensor arrays deployed in water [154] or in ice [155] rely on the acoustic 
detection technique [156, 157]. 

The radio detection technique is traditional at the South Pole, the RICE detector [147] 
operating for a long time has paved the path toward the ARA project [149] consisting 
of a series of antennas in the firn ice layer. The ARA sensitivity to ultra high energy 
neutrinos would eventually reach the level of unfavorable GZK neutrino flux expectations 
(pure Fe composition). ARA could be complemented with surface radio stations [158] 
and, conceivably with acoustic detection devices. 

Contrary to salt and water, ice has the unique feature of allowing the detection of three 
concurrent signatures accompanying a neutrino event, optical Cherenkov light [159], co- 
herent radio Cherenkov and thermoacoustic emissions. The combined detection of two or 
more of these signatures, leading to a strong background reduction, would unambiguously 
tag a neutrino event and thus firmly establish the event origin [160, 161]. 

In ice, the attenuation length of the radio and acoustic emission are respectively 
O(km) and 0(0.3 km) [162, 163, 147, 164]. The applicability of the acoustic detection 
technique from a physical standpoint is however not so much related to the attenuation 
length: an inadequately high absolute noise level may definitely impede its applicability 
in an inappropriate environment. This is still an unknown issue and solutions to remedy 
have been proposed [165]. 



Acoustic detection method. The emission of a thermo- 
acoustic pulse following a UHE neutrino interaction is at 
the origin of the detection technique: the sudden energy 
release in a medium from the particles along the shower 
provokes its thermal expansion with the subsequent shock 
wave emission. The acoustic pulse is coherent normal to the 
shower axis and therefore exhibits a "pancake- like" shape 
as illustrated Fig. 16. It can be described by means of the 
wave equation that follows from the equations of motion, 
continuity and state, 



(63) 



V 2 (p- 



1 1 .. 

—P) -~2P = 

w c z 



P_dE_ 
C~ p ~dt 




which includes a source term on the r.h.s. and where p is 

the pressure, ujq = ^o(f) a characteristic attenuation fre- 16- - UHE neutrino interac- 



quency, which depends on medium compressibility and vis 
cosity and frequency (often neglected), c the sound speed, (3 
the bulk coefficient of thermal expansion and C p the specfic 
heat at constant p. The heat deposition is quasi instanta- 
neous, i.e. pe{*' ,t) = pE{r')S{t) reducing the problem to 
an integral solution, 



tion radio, acoustic and optical 
signatures. 



(64) 



p(r,t)= J p E (r')G(r-r',t)dr' 
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with 

, . , 8 t-r/c (t-r-/cf 

where t = y/r/ {loqc). G(r, t) is proportional to the pressure pulse generated by a point- 
like energy deposition Pe{ v ) = Po^ ( r ) with an apparent bipolar shape. 

In the typical media (salt, ice, water), the peak frequencies of the acoustic pulse are 
in the range 10-50 kHz, i.e. the signal remains coherent (most of the energy deposition 
is within a distance of order of a few cm from the shower axis, corresponding to typical 
frequencies up to 0(100 kHz). 

The efficiency of pressure wave generation is described by the Griineisen medium 
parameter j g = Be 2 /C p . The c factor appears when integrating expression (64) for a 
shower- like energy distribution pe (r) and translates the higher frequencies of the coherent 
emission, which can be attained. j g is about 7 times larger in ice than in seawater due 
to the sound speed difference Ci cc /c wa t or ~ 2.5. Ratios of bulk coefficient of thermal 
expansion to specific heat are roughly equal for water and ice( 6 ). Coincidentally, in ice, 
100 kHz is the predicted value above which the ice attenuation length due to scattering 
at crystal boundaries quickly drops [162]. The acoustic detection technique in general 
and its application in seawater is discussed in detail in [166]. 

Experimentally, the layout of an acoustic array would would be composed of sparse 
strings densely (pancake shape) instrumented. In the ice, it would be deployed down 
to depth smaller than about 1 km, as the absorption length acutely depends on the ice 
temperature (and the temperature gradient is due to the geothermic flux). 



Radio detection method. The charge asymmetry in the development of an EM shower 
and the subsequent coherent radio emission (a short linearly polarized pulse of about a 
ns) is at the origin of the detection technique. A shower develops an approx. 20% charge 
a symmetry due roughly to the combined effect of shorter pathlcngths of the positrons 
(in-flight anihilation) and of an amplification of the electronic component, via ionization 
processes (Compton, M0ller, Bhabha), 

1. Shower photons interact with atomic electrons via Compton scattering, the main 
mechanism for the generation of the charge assymetry, 7 + e~ tom — > 7 + e~ 

2. Shower positrons interact with atomic electrons via Bhabha scattering, e + e~ tom — > 
e + e~. This process inject new energetic electrons in the shower while reducing the 
incident positron energy. 

3. Via anihilation in-flight, which terminates an incident positron trajectory, e + e~ tom — > 

77 

4. Shower electrons interact with atomic electrons via M0ller scattering): e~e~ tom — > 

c~c~ 

All the rclativistic particles in the shower contribute to the radio Cherenkov emission 
(of order of 10 14 particles in a 10 EeV shower). While the Cherenkov emission is faint 



( 6 ) c icc = 3.92 • 10 3 [m/s], c watcr = 1.53 • 10 3 [m/s], f3 icc = 1.25 • 10" 4 m 3 m" 3 K"\ /? wator = 
2.55- 10" 4 m 3 m- 3 K- 1 C p ,ico = 1.72- 10 3 J kg" 1 KT 1 , C p , wat er = 3.9- 10 s Jkg^K -1 , see [163]. 
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at radio wavelength, the amplification of the signal is due to the coherence of the col- 
lective emission. The coherence compensates the suppression of this emission at long 
wavelengths and for L of order 10 cm shower core lateral extension. Coherence up to 
about A > 4L is expected, i.e. up to about GHz frequencies. 

Experimentally, the Cherenkov radio frequency signal is exploited for UHE neutrino 
detection by operating antennas in the range 0.1 - 1 GHz. The South Pole ice has 
adequate properties for the propagation of radio waves at depths down to 1.5 km with 
an absorption length of 1.2 km (at a few 100's of Mhz). 

12. Conclusions 

From instrumental, phenomenological and methodological standpoints, neutrino as- 
tronomy has come a long way since the 60's. In these lectures, we glimpsed at the 
methodology in a broad sense, trying to establish a link to phenomenology, present and 
future instrumental techniques. 

While the origin of the CRs remains a main focus for neutrino telescopes, the research 
scope has widened to include questions related to fundamental particle physics and cos- 
mology. Neutrino telescopes offer a rich and promising research program for the next 
decades. 

* * * 

I thank the students and organizers and fellow lecturers for a nice time at the school. 
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